Modeling Urban Growth Patterns with Correlated Percolation 
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We propose and test a model that describes the morphology of cities, the scaling of the ur- 
ban perimeter of individual cities, and the area distribution of systems of cities. The model is 
also consistent with observable urban growth dynamics, our results agreeing both qualitatively and 
quantitatively with urban data. The resulting growth morphology can be understood from interac- 
tions among the constituent units forming an urban region, and can be modeled using a correlated 
percolation model in the presence of a gradient. 



I. INTRODUCTION 

Traditional approaches to urban science as exempli- 
fied in the work of Christaller jlj, Zipf 0h Stewart and 
Warntz || , Beckmann [Q , and Krugman ||| are based on 
the assumption that cities grow homogeneously in a man- 
ner that suggests that their morphology can be described 
using conventional Euclidean geometry. However, recent 
studies have proposed Q that the complex spatial phe- 
nomena associated with actual urban systems is rather 
better described using fractal geometry consistent with 
growth dynamics in disordered media [0-^). 

Predicting urban growth dynamics also presents a chal- 
lenge to theoretical frameworks for cluster dynamics in 
that different mechanisms clearly drive urban growth 
from those which have been embodied in existing phys- 
ical models. In this paper, we develop a mathematical 
model that relates the physical form of a city and the 
system within which it exists, to the locational decisions 
of its population, thus illustrating how paradigms from 
physical and chemical science can help explain a uniquely 
different set of natural phenomena - the physical arrange- 
ment, configuration, and size distribution of towns and 
cities. Specifically, we argue that the basic ideas of perco- 
lation theory when modified to include the fact that the 
elements forming clusters are not statistically indepen- 
dent of one another but are correlated, can give rise to 
morphologies that bear both qualitative and quantitative 
resemblance to the form of individual cities and systems 
of cities. Some of these results are briefly described in 
Ref. @. 

We consider the application of statistical physics to ur- 
ban growth phenomena to be extremely promising, yield- 
ing a variety of valuable information concerning the way 
cities grow and change, and more importantly, the way 
they might be planned and managed. Such information 
includes (but is not limited to) the following: 

(i) the size distributions of towns, in terms of their 
populations and areas; 

(ii) the factal dimensions associated with individual 



cities and entire systems of cities; 

(Hi) interactions or correlations between cities which 
provide insights into their interdependence; 

(iv) the relevance and effectiveness of local planning 
policies, particularly those which aim to manage and con- 
tain growth. 

The size distribution of cities has been a fundamental 
question in the theory of urban location since its incep- 
tion in the late 19th century. In the introduction to his 
pioneering book published over 60 years ago, Christaller 
|l| posed a key question: "Are there laws which deter- 
mine the number, size, and distribution of towns?" This 
question has not been properly answered since the pub- 
lication of Christaller's book, notwithstanding the fact 
that Christaller's theory of central places [0 and its elab- 
oration through theories such the rank-size rule for cities 
[Q-Q embody one of the cornerstones of human geogra- 
phy. 

Our approach produces scaling laws that quantify such 
distributions. These laws arise naturally from our model, 
and they are consistent with the observed morphologies 
of individual cities and systems of cities which can be 
characterized by a number of fractal dimensions and per- 
colation exponents. In turn, these dimensions are con- 
sistent with the density of location around the core of 
any city, and thus the theory we propose succeeds in ty- 
ing together both intra- and inter-urban location theories 
which have developed in parallel over the last 50 years. 
Furthermore, the striking fact that cities develop a power 
law distribution without the tuning of any external pa- 
rameter might be associated with the ability of systems 
of cities to "self-organize" [|| . 



II. THE DLA MODEL 



Cities grow in a way that might be expected to resem- 
ble the growth of two-dimensional aggreg ates of particles, 
and this has led recent attempts ||,[ll],|l2| to model urban 
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growth using ideas from the statistical physics of clus- 
ters. In particular, the model of diffusion limited aggre- 
gation (DLA) Jl|[14[ has been applied to describe urban 
growth [pj, and results in tree- like dendritic structures 
which have a core or "central business district" (CBD). 
The DLA model is a physical model used to describe ag- 
gregation phenomena and is related to problems from the 
field of oil recovery in which "viscous fingering" occurs 
when a low viscosity fluid is pushed under pressure into a 
fluid with a larger viscosity (as occurs when an oil field is 
flooded with water in an attempt to "push out the oil"). 

The DLA model predicts that there exists only one 
large fractal cluster that is almost perfectly screened 
from incoming "development units" (people, capital, re- 
sources, etc), so that almost all the cluster growth oc- 
curs in the extreme peripheral tips. However, quantita- 
tive data do not support all the properties of the DLA 
model. For instance, the DLA model predicts that the 
urban population density p(r) decreases from the city 
center as a power law, 



p(r) 



„D-2 



(1) 



where r is the radial distance from the core, and D ~ 1.7 
is the fractal dimension of DLA. However, urban data 
have been more commonly fit to an exponential decay 
]l5fl . In the DLA model only one large central place or 
cluster is generated, while a real urban area is formed by 
a system of central places that are spatially distributed in 
a hierarchy of cities. Still another concern regarding the 
morphology of the DLA model is that DLA is a simply- 
connected cluster. Cities grow in a more compact way, 
with a well-defined urban boundary or external perime- 
ter not accounted for by the dendritic fractal growth of 
DLA. 

Here we show that an alternative model, in which de- 
velopment units are correlated rather than being added 
to the cluster at random, is better able to reproduce the 
observed morphology of cities and the area distribution 
of sub-clusters ("towns") in an urban system, and can 
also describe urban growth dynamics. Our "physical" 
model 10 1, which corresponds to the correlated perco- 
lation model |l6]-|o| in the presence of a density gradi- 
ent [pT|~^3|, is motivated by the fact that in urban areas 
development attracts further development. The model 
offers the possibility of predicting the global properties 
(such as scaling behavior) of urban morphologies. 



III. CORRELATED PERCOLATION MODEL 

In the model we now develop, we take into account two 
points: 

(i) First, data on population density p(r) of actual ur- 
ban systems are known to conform to the relation |15| 



p(r) = Pee 



-\r 



(2) 



where r is the radial distance from the central business 
district (CBD) situated at the core, and A is the density 
gradient. The density gradient quantifies the extent of 
the urban spread around the central core. The probabil- 
ity that a unit occupies a given spot decreases gradually 
as the distance from a central, compact core increases. 

(ii) Second, in actual urban systems, the development 
units are not positioned at random. Rather, there exist 
correlations arising from the fact that when a develop- 
ment unit is located in a given place, the probability of 
adjacent development units increases naturally; each site 
is not independently occupied by a development unit, 
but is occupied with a probability that depends on the 
occupancy of the neighborhood. In urban settings, devel- 
opment units do not attach themselves randomly to an 
existing cluster. Their placement is strongly influenced by 
the presence of other units. When a unit occupies a cer- 
tain location, the probability of additional development 
is highest in its vicinity, and this probability decreases 
at a certain rate as the distance from the unit increases. 
Thus, the rules of placement are affected by long-range 
"interactions" that influence how clusters form and grow. 
What happens at a given site depends on the state of 
many other sites. These correlations reflect the tendency 
of people to locate next to one another, as articulated in 
traditional urban science as economies of urban agglom- 
eration. 

In order to quantify these ideas, we consider the cor- 
related percolation model |16|- P0| in the presence of a 
concentration gradient |2l|-p3|. We start by describing 
the uncorrelated site percolation problem, which corre- 
sponds to the limit where correlations are so small as 
to be negligible |7|-[|. We first define a random number 
u{r), called the occupancy variable, at every site r=(i,j) 
in a square lattice of L 2 sites. The numbers u(r) are un- 
correlated numbers with a uniform probability distribu- 
tion between [0, 1]. A site in the lattice is occupied if the 
occupancy variable u(r) is smaller than the occupation 
probability p, which is a quantity fixed for every site in 
the lattice. A cluster is a set of sites connected via near- 
est neighbor sites. When p is small only isolated clusters 
exist. At a critical value of the concentration called p c 
an "incipient infinite cluster" forms which, for a finite 
system, connects two sides of the system. 

Our basic model is a percolation model modified to in- 
troduce correlations among the units and the fact that 
the concentration p is not constant for all the points in 
the lattice but presents the gradient shown in Eq. (Q) . In 
our model we consider "development units" represent- 
ing buildings, people, resources which are added to the 
cluster in similar fashion as in percolation. Since these 
units are added in a correlated fashion, we next consider 
a modification of the percolation problem to incorporate 
correlations among the occupancy variables u(r). 

To introduce correlations among the variables we use 
a method proposed in |2(J which is a modification of 
the Fourier filtering method (Ffm) p^-p6 18 suitable 
for large systems. Consider a stationary sequence of L 2 
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uncorrelated random numbers {u(r)}, r = = 
1,..,L. The correlation function is (u(r) u(r')) ~ 5 r y , 
with 5 r y the Kronecker delta, and the brackets denote an 
average with respect to a Gaussian distribution. We use 
the sequence {u(r)}, to generate a new sequence, {r](r)}, 
with a long-range power-law correlation function C (£) of 
the form BO) 



C(t) = (u(r) u{r')) = {l+l 2 )- a ' 2 , 



(3) 



where, i — \r — r'\, a is the correlation exponent, and the 
long-range correlations are relevant for < a < d = 2, 
where d is the dimension of the substrate — a > 2 cor- 
responds to the uncorrelated problem, and a — > to the 
strongly correlated problem. 

The spectral density S(q), defined as the Fourier trans- 
form of C(£), has the form 



S(q) 



2n 



i) 



q\02 
2, 



(4) 



where q — \q\, qi — 2itmi/ L, —L/2 < rtii < L/2, i = 1,2, 
and P2 = {ot — 2)/2. {r)(q)} are the Fourier transform 
coefficients of {r)(r)}, and satisfy 



V (q) = (S(q)) 1/2 u(q) 



(5) 



where {u(q)} are the Fourier transform coefficients of 
Mr)}. 

The actual numerical algorithm for Ffm consists of the 
following steps: (i) Generate a two-dimensional sequence 
{u(r)} of uncorrelated random numbers with a Gaussian 
distribution, and calculate the Fourier transform coeffi- 
cients {u(q)}. (ii) Obtain {r/(q)} using (||) and (||). (Hi) 
Calculate the inverse Fourier transform of {i](q)} to ob- 
tain {rj(r)}, the sequence in real space with the desired 
power-law correlation function which asymptotically be- 
haves as 



c(i) - r 



(G) 



The assumption of power-law interactions is motivated 
by the fact that the "decision" for a development unit to 
be placed in a given location decays gradually with the 
distance from an occupied neighborhood. 

Finally we consider that the development units are po- 
sitioned with a probability which behaves in the same 
fashion as known for cities Eq. (||). Therefore we relax 
the assumption that the concentration p is constant for 
all the points in the lattice, and we consider that the de- 
velopment units are positioned with an occupancy prob- 
ability 



p(r) = p(r)/p , 



(7) 



that behaves in the same fashion as is known in observa- 
tions of real cities. This last modification corresponds to 
the percolation problem in the presence of a concentra- 
tion gradient proposed in |2l|-|25|]. 



A. 





FIG. 1. Simulations of urban systems for different degrees 
of correlation. Here, the urban areas are red, and the exter- 
nal perimeter or urban boundary of the largest cluster con- 
nected to the CBD is light green. In all the figures, we fix 
the value of the density gradient to be A = 0.009. (a) and 
(b) Two different examples of interactive systems of cities for 
correlation exponents a — 0.6 and a = 1.4, respectively. The 
development units are positioned with a probability that de- 
cays exponentially with the distance from the core. The units 
are located not randomly as in percolation, but rather in a 
correlated fashion depending of the neighboring occupied ar- 
eas. The correlations are parametrized by the exponent a. 
The strongly correlated case corresponds to small a (a — > 0). 
When a > d, where d is the spatial dimension of the sub- 
strate lattice (d = 2 in our case) , we recover the uncorrelated 
case. Notice the tendency to more compact clusters as we in- 
crease the degree of correlations (a — > 0). (c) As a zeroth 
order approximation, one might imagine the morphology pre- 
dicted in the extreme limit whereby development units are 
positioned at random, rather than in the correlated way of 
Figs, la and 16. The results for this crude approximation of a 
non-interactive (uncorrelated) system of cities clearly display 
a drastically different morphology than found from data on 
real cities (such as shown in Fig. 0a,) . The non-interactive limit 
looks unrealistic in comparison with real cities, for the lack of 
interactions creates a urban area characterized by many small 
towns spread loosely around the core. 
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In order to apply the above procedure to the percola- 
tion problem, we study the probability distribution P(rj) 
of the correlated sequence rj(r). We find that when the 
uncorrelated variables u(r) are taken from a Gaussian 
distribution, rj{r) also has a Gaussian distribution. We 
next discretize the variables generating a sequence , 
according to = Q(6—rj(r)) where 9 is chosen to sat- 

isfy p(r) = j_ P(rj)dr), with p{r) the occupancy proba- 
bility and is the Heaviside step function. 

Notice that we have defined two different properties 
of the system. First we introduced long range correla- 
tions among the variables by modifying the occupancy 
variables r)(r). These correlations are isotropic, i.e., all 
the points in space are connected by interactions quanti- 
fied by a power law. The fact that we consider a slowly- 
decaying power-law scale-free function is due to the fact 
that any other correlation function will display a cut-off 
after which correlations are negligible. Since we are look- 
ing at properties of actual cities at large length scales, a 
coarse grain will transform a finite range correlated sys- 
tem into a uncorrelated system, i.e., a system with a finite 
cut-off in the correlations becomes uncorrelated at large 
scales. This situation occurs when we consider power law 
correlations of the form (Q) since it is a scale-free func- 
tion. Thus correlations are expected to be relevant at all 
length scales. One must distinguish the type of correla- 
tion introduced by (||) from the correlations arising at 
the critical concentration p c . In this case, the connected- 
ness length of the system is said to be infinite since two 
occupied sites separated by an arbitrary distance may be 
connected by the infinite cluster, and thus correlated in 
space. However, the correlations introduced by (||) goes 
beyond this type of connection between sites. Due to cor- 
relation of type (||), even sites which belong to different 
clusters are correlated. 

Second, we consider that the probability of occupancy 
of the sites decays exponentially with the center point 
always occupied. This property of the system is indepen- 
dent of the type of correlation chosen. The correlation 
exponent a and the density gradient A are the only pa- 
rameters of the model to be determined by empirical ob- 
servations. 



IV. STATICS 

We first discuss the influence of the correlations on the 
morphology of a system of cities generated in the present 
model. Therefore we fix the value of the concentration 
gradient A in Eq. (^) and we show in Fig. |l| our simula- 
tions of correlated urban systems for different degrees of 
correlation. We see that the larger the degree of correla- 
tions the more compact the clusters are. The correlations 
have the effect of agglomerating the units around an ur- 
ban area. In the simulated systems the largest city is sit- 
uated in the core (which acts as the "attractive" center 
of the city), and this is surrounded by small clusters or 



"towns." The correlated clusters are fairly compact near 
their respective centers and become less compact near 
their boundaries, in qualitative agreement with empiri- 
cal data on actual large cities such as Berlin, Paris and 
London, (see i.e. Refs. pj,^). The strongly correlated 
case of Fig. |l]a a — > results in a system of cities look- 
ing more realistic than the uncorrelated case (Fig. |l|c). 
The uncorrelated case results in a systems of very small 
cities spread around a central city, while the cities in the 
correlated case look more compact and more realistic. 




FIG. 2. Qualitative comparison between the actual urban 
data and the proposed model, (a) Three steps of the growth 
with time of Berlin and surrounding towns. Data are shown 
for the years 1875, 1920, and 1945 (from top to bottom), (b) 
Dynamical urban simulations of the proposed model. We fix 
the value of the correlation exponent to be a — 0.05 (strongly 
correlated case), and choose the occupancy probability p(r) 
to correspond to the density profiles shown in Fig. ^[ We use 
the same seed for the random number generator in all figures. 

For distances smaller that rt, there is a high concen- 
tration of sites since p(r) > p c , and the cluster thus gen- 
erated plays the role of the infinite cluster. For distance 
larger than r/, we have p(r) < p c , so that only isolated 
clusters exist, which form the system of small cities sur- 
rounding the large city situated in the core. 

The urban boundary of the largest city is defined to 
be the external perimeter of the cluster connected to 
the CBD. Since p(r) decreases as one moves away from 
the core, the probability that the largest cluster remains 
connected decreases with r. The mean distance of the 
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perimeter from the center rf is determined by the value 
of r for which p(r) equals the percolation threshold — i.e. 
p(r f ) =p c , so 0-61 



L(£) 



l>D c 



(12) 



Tf = A -1 ln(l/p c ). 



(8) 



The geometrical properties of the external perimeter 
of the largest city correspond to the properties of the ex- 
ternal perimeter of the infinite cluster of the percolation 
problem in the absence of a gradient [|T). The critical 
properties of the clusters can be analyzed in terms of the 
percolation properties. Percolation clusters formed be- 
low p c are characterized by a finite connectedness length 
which is the typical distance at which two sites are ex- 
pected to be connected via nearest neighbor sites (do 
not confuse with the correlations introduced via Eq. (ph). 
This connectedness length diverges when the infinite clus- 
ter forms at p c , i.e. £ ~ \p — p c \ ~ v , where v is the con- 
nectedness length exponent. In the case of gradient per- 
colation the clusters formed below p c for r > rf are char- 
acterized by this length, which is now a function of the 
distance r 



£(r) - \p(r) -p c \" 



(9) 



Moreover, due to the existence of long range correlation 
among the variables the exponent v is not universal, but 
changes continuously with the degree of correlation given 
by a We will see that several critical properties of 
the percolation clusters change with the correlations. 

The width Of of the external perimeter of the largest 
city is defined as 



a f ^((r-r f )Y /2 , 



(10) 



where r/ = (r), and r belongs to the external perimeter 
of the central cluster. The width of the external perime- 
ter is a function of the concentration gradient A and it is 
known to scale as [Ell 



(11) 



The value of v corresponding to the uncorrelated per- 
colation problem is v = 4/3. However the presence of 
long range correlations of the type (0) drastically affects 
the value of the connectedness exponent, which is now a 
function of a, v(a) as observed in previous studies of long 
range correlated percolation [Jl7|,[l8|. We have simulated 
the correlated percolation problem with a gradient and 
using Eq. (|ll| ) we find a drastic increase of v(a) with the 
increase of the long-range correlations (a — > 0) (Figs. |^a, 
|§|b). In particular v(a) seems to increase very drastically 
for a system of strong correlations a — > 0. In fact for such 
a system, we expect a mean field situation where all sites 
in the lattice are connected to the rest of the sites. In 
this case the percolation threshold for the site percola- 
tion problem should be p c = 0.5 and the connectedness 
length should be zero below p c and infinite above p c . 

The scaling of the length of the urban boundary of the 
largest city within a region of size I, 



defines the fractal dimension D e , which we calculate to 
have values D e ~ 1.33 for the uncorrelated case, and 
D e ~ 1.4 for strong correlations (a — > 0) (Fig. ||c). The 
small variation of the fractal dimension of the external 
perimeter does not rule out the fact that it may be inde- 
pendent of the correlations. These values are consistent 
with actual urban data, for which values of D e between 
1.2 and 1.4 are measured ||. 

Near the frontier and on length scales smaller than the 
width of the frontier a/, the largest cluster has fractal 
dimension df ~ 1.89, as defined by the "mass-radius" 
relation 



M(r) 



(13) 



where M{r) is the mass of the cluster inside a region of 
radius r. The value df — 1.89 corresponds to the fractal 
dimension of the uncorrelated percolation clusters and 
we find that it is valid independent of the correlations. 
However, as a — > we expect a compact cluster with 
dimension df = 2. The fact that we are unable to see 
this limit might be due to numerical limitations near the 
mean field point a = 0. 

The number of sites of the frontier Nf also scales with 
the concentration gradient E| 



N 



f 



(14) 



This relation provides another way of calculating the 
fractal dimension df and the exponent v, which we used 
to verify our calculations. 

It is important to stress that under the present perco- 
lation picture cities are fractal structures only near the 
external perimeter of the largest city, and on length scales 
smaller than the width of the frontier defined by Eq. (|Tc|). 
The width is a function of the concentration gradient A, 
Eq. (11) so that the larger the spread of the city the 
larger the region where the city is fractal. However, at 
distances close to the center of the largest city, the clus- 
ter is clearly non fractal since p(r) > p c and the cluster 
becomes compact. On the other side for larger distances 
p(r) < p c , only small isolated clusters exist with a def- 
inite connectedness length associated with them (||), so 
that they are not fractal either. 



V. AREA DISTRIBUTION OF URBAN 
SETTLEMENTS 

So far, we have argued how correlations between oc- 
cupancy probabilities can account for the irregular mor- 
phology of towns in a urban system. As can be seen in 
Fig. Ep, the towns surrounding a large city like Berlin 
are characterized by a wide range of sizes. We are inter- 
ested in the laws that quantify the town size distribution 
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N(A), where A is the area occupied by a given town or 
"mass" of the agglomeration. 



X. 10' 

O 



^x», 

• a=0.2 
□ a=0.8 

♦ 01=1.4 
a uncorrelated 



(a) 



X 



and Berlin for different years (Fig. |]a), and also we an- 
alyzed the area distribution of urban systems at larger 
scales by using the data of all settlements of Great Britain 
for the years 1981 and 1991 |28|| . In the case of the towns 
around Berlin and London, we first digitize the empirical 
data of Fig. 4.1 of Ref. |f| (Berlin 1920 and 1945, shown 
in the last two panels of Fig. and Fig. 10.8 of Ref. 
H corresponding to London 1981. Then, we count the 
number of towns that are covered by A sites, putting the 
result in logarithmically spaced bins (of size 1.2 fe , with 
k = 1, 2, ...), and averaging over the size of the bin. 

We calculate the actual distribution of the areas of the 
urban settlements around Berlin and London, and find 
(Fig. ^p.) that for both cities, N(A) follows a power-law 
with exponent close to 2 



N{A) ~ A 



-1.98 



[Berlin, 1920, 1945], 



(15) 




N(A) ~ A 



-1.96 



[London, 1981]. 



(16) 



FIG. 3. (a) Width of the external perimeter as a function 
of the density gradient, 07(A), for several degree of correla- 
tions, (b) Connectedness length exponent v(a) as a function 
of the correlation exponent a calculated from Fig. ^ using 
Eq. ([n]). The value u = 4/3 corresponds to the uncorrelated 
percolation problem (a = 2). (c) Fractal dimension of the 
external perimeter of the largest cluster as a function of the 
degree of correlation, D e (a). 



We have analyzed the distribution of areas of actual 
cities, such as the system of cities surrounding London 



Figure ^3 shows the distribution of all urban areas in 
Great Britain for the years 1981 and 1991.. We find a 
power-law with an exponent consistent with the data of 
London and Berlin at smaller scales: 



N(A) ~ A 



-2.03 



[Great Britain, 1981, 1991]. (17) 



Other studies have recently confirmed the validity of 
these results for larger length scales and also for the pop- 
ulation distributions which is known to scale as the oc- 
cupied area p9fl . 

These results can be understood in the context of our 
model. Insight into this distribution can be developed 
by first noting that the small clusters surrounding the 
largest cluster are all situated at distances r from the 
CBD such that p(r) < p c or r > r/. Therefore, we find 
N(A), the cumulative area distribution of clusters of area 
A, to be 



N(A) 



Here. 



n(A,p) dp- A-( T+1 ' d f v \ (18) 



n(A,p)~A- r g(A/A ) (19) 



is defined to be the average number of clusters con- 
taining A sites for a given p at a fixed distance r, and 
t = 1 + 2/df. Here, 



Mr)~Z(r) d f~\p(r)-p c \- d f 



(20) 



corresponds to the maximum typical area occupied by 
a cluster situated at a distance r from the CBD, while 
g(A/Ao) is a scaling function that decays rapidly (expo- 
nentially) for A > Aq. 

In our numerical simulations we find a drastic increase 
of v{a) with the increase of the long-range correlations 
(a — > 0) (Fig. ||b) The connectedness exponent v{a) af- 
fects the area distribution of the small clusters around the 



G 



CBD (Fig. @), as specified by Eq. @, and can be used to 
quantify the degree of interaction between the CBD and 
the small surrounding towns. For instance, for a strongly 
correlated system of cities characterized by small a, v{a) 
is large so that the area A${r) and the linear extension 
£(r) of the towns will be large even for towns situated 
away from the CBD. This effect is observed in the corre- 
lated systems of cities of Fig. [j]. 

In Fig. |5^ we plot the power-law for the area distribu- 
tion predicted by Eq. ( |hS| ) along with the real data for 
Berlin and London and all Great Britain. In particular, 
the slope predicted for the uncorrelated system is 



(21) 



N(A)~A 2 45 , [uncorrelated model], 

while for the strongly correlated model is 

N(A) ~ A~ 2 - 06 , [strongly correlated model]. (22) 

Therefore, we find that the power law of the area distri- 
bution of actual cities are consistent with the prediction 
(dashed line, Fig. [|a) for the case of highly correlated 
systems. These results quantify the qualitative agreement 
between the morphology of actual urban areas and the 
strongly correlated urban systems obtained in our sim- 
ulations. Clearly, the exponent of the area distribution 
provides a stronger test of our model against observations 
than does the fractal dimension D e of the perimeter. 
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remarkable pattern of decentralization, which is quanti- 
fied by the decrease of A(t) with time (see Table 4 in Ref. 
p0|] , and Fig. |J). Therefore the dynamics in the model 
are quantified by a decreasing A(t), as occurs in actual 
urban areas. In the context of our model, this flattening 
pattern can be explained as follows. The model of perco- 
lation in a gradient can be related to a dynamical model 
of units (analogous to the development units in actual 
cities) diffusing from a central seed or core ]2l|-|23|] . In 
this dynamical system, the units are allowed to diffuse on 
a two-dimensional lattice by hopping to nearest-neighbor 
positions. The density of units at the core remains con- 
stant: whenever a unit diffuses away from the core, it is 
replaced by a new unit. 
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FIG. 4. Log-log plot of the area distribution function N(A) 
calculated for the present model for different degrees of cor- 
relation. From top to bottom, a = 0.2, a = 0.8, a = 1.4, and 
uncorrelated case. The linear fits correspond to thepredic- 
tions of Eq. ( |l8| ) using the values of u(a) from Fig. |^b, and 
d f = 1.89. 



VI. DYNAMICS 



We now discuss a generalization of our static model to 
describe the dynamics of urban growth. Empirical stud- 
ies Q of the population density profile of cities show a 



7 



FIG. 5. (a) Log-log plot of the area distribution N(A) of 
the actual towns around Berlin and London. We first digi- 
tize the empirical data of Fig. 4.1 of Ref. @ (Berlin 1920 
and 1945, shown in the last two panels of Fig. Bp), and Fig. 
10.8 of Ref. U (London 1981). Then, we count the number of 
towns that are covered by A sites, putting the result in log- 
arithmically spaced bins (of size 1.2 fc , with k = 1,2, ...), and 
averaging over the size of the bin. A power-law is observed 
for the area distributions of both urban systems. The dotted 
line shows the predictions of our model for the uncorrelated 
case (slopes 2.45), while the dashed line gives results for the 
strongly correlated case (slope= 2.06). Note that the area dis- 
tributions for both cities agree much better with the strongly 
correlated case (a — > 0). (b) Log-log plot of the area distribu- 
tion of all the urban areas in Great Britain in 1981 and 1991. 
The data fits to a power law of exponent 2.03. Notice also the 
very small changes of the urban areas in a 10 year period. 

The density of units can be mapped to the density of 
occupied urban areas 



Pa{t) 



-Xr 



(23) 



which in turn is proportional to the population density 
p(r) ||. A well-defined diffusion front, defined as the 
boundary of the cluster of units that is linked to the cen- 
tral core, evolves in time. The diffusion front corresponds 
to the urban boundary of the central city. The static 
properties of the diffusion front of this system were found 
to be the same as those predicted by the gradient per- 
colation model plj-p^]. Moreover, the dynamical model 
can explain the decrease of X(t) with time observed em- 
pirically. As the diffusion front situated around rt moves 
away from the core, the city grows and the density gra- 
dient decreases since X(t) oc 1/V/. 

These considerations are tested in Fig. which shows 
our dynamical urban simulations of a strongly interacting 
system of cities characterized by a correlation exponent 
a = 0.05 for three different values of A obtained from 
the data of Berlin from Fig. || Qualitative agreement is 
observed between the morphology of the cities and towns 
of the actual data of Fig. Wn and the simulations of Fig. 
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FIG. 6. Semi-log plot of the density of occupied urban ar- 
eas Pa{t) — e~ Ar for the three different stages in the growth 
of Berlin shown in Fig. ^a. Least square fits yield the results 
A ~ 0.030, A ~ 0.012, and A ~ 0.009, respectively, showing 
the decrease of A with time. We use these density profiles in 
the dynamical simulations of Fig. 



VII. DISCUSSION: URBAN PLANNING 

Throughout this century, the dominant planning pol- 
icy in many western nations has been the containment 
of urban growth. This has been effected using several 
instruments, particularly through the siting of new set- 
tlements or new towns at locations around the growing 
city which are considered to be beyond commuting dis- 
tance, but also through the imposition of local controls 
on urban growth, often coordinated regionally as "green 
belts" |3l| . One of the key elements in the growth models 
we have proposed here is the characteristic length scale 
over which growth takes place. In the case of the gradi- 
ent percolation model, correlations occur over all length 
scales, and the resulting distributions are fractal, at least 
up to the percolation threshold. 

In examining the changing development of Berlin in 
Fig. ||a, it appears that the fractal distribution remains 
quite stable over a period of 70 years and this implies that 
any controls on growth that there might have been do 
not show up in terms of the changing settlement pattern, 
implying that the growth dynamics of the city are not 
influenced by such control. A rather different test of such 
policies is provided in the case of London where a green 
belt policy was first established in the 1930s and rigor- 
ously enforced since the 1950s. The question is whether 
this has been effective in changing the form of the set- 
tlement pattern. First, it is not clear that the siting of 
new towns beyond London's commuting held was ever 
beyond the percolation field and thus it is entirely pos- 
sible that the planned new settlements in the 1950s and 
1960s based on existing village and town cores simply 
reinforced the existing fractal pattern. 

In the same manner, the imposition of local controls 
on growth in terms of preserving green field land from 
development seems to have been based on reinforcing 
the kind of spatial disorder consistent with morphologies 
generated through correlated percolation. The regional 
green belt policy was based on policies being defined lo- 
cally and then aggregated into the green belt itself, and 
this seems to suggest that the morphology of nondevel- 
opment that resulted was fractal. This is borne out in 
a fractal analysis of development in the London region 
which suggests that the policy has little impact on the 
overall morphology of the area ||[32]]. Moreover, we note 
that the coincidence between the settlement area distri- 
bution for different cities and different years (Berlin 1920 
and 1945, and London 1981) suggests that local planning 
policies such as the green belt may have a relatively low 
impact on the distribution of towns. Our model suggests 
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that the area distribution is determined by the degree of 
interactions among development units, and that its scal- 
ing properties are independent of time. Current debates 
on urban growth have now shifted to the development of 
brownfield sites in cities, and it would be interesting to 
quantify the extent to which such future developments 
might reinforce or counter the "natural" growth of the 
city as implied in these kinds of models. 

To develop more detailed and conclusive insights into 
the impact of urban policies on growth, it is necessary 
to develop the model further. This model implies that 
the area and size distributions, the degree of interaction 
amongst dependent units of development, and fractal di- 
mension are independent of time. The only time depen- 
dent parameter is the gradient A and it appears that we 
might predict future urban forms simply by extrapolat- 
ing the value of A in time. However, we have yet to in- 
vestigate the influence of topography and other physical 
constraints on development, the influence of transport 
routes and the presence of several "independent" central 
cores or CBDs in the urban region. 

These models can also be further adapted to predict 
bond as well as site percolation and in future work we 
will explore the extent to which such interactions be- 
tween sites and cities might be modeled explicitly. Our 
interest in such examples is in the universality of the ex- 
ponents that we have demonstrated here, and which we 
wish to relate to the impact of urban planning policies. 
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